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I. INTRODUCTION 

In this paper, we develop an exact analytical 
representation of neutrino oscillations jl]] in mat¬ 
ter within the framework of the Standard Neutrino 
Model (SNM) [2| with 3 Dirac Neutrinos. The ex¬ 
act closed-form expressions we give for the time- 
evolution operator S(t,t') are obtained from H u 
using the Lagrange interpolation formula given in 
Ref. Q- The resulting expressions are easily evalu¬ 
ated without any approximations. 

The paper is divided into two main parts. In 
the first, we summarize the main results of our the¬ 
ory. Details underlying the derivation are given in 
Appendices. We also retrieve the well-known two- 
neutrino flavor results as a special case of our general 
results. 

In the second part we address other analytical 
formulations found in the literature. The expan¬ 
sion of the neutrino oscillation probability in one of 
the small parameters a and sin 2 $13 of the standard 
neutrino model (SNM) for H u is of particular inter¬ 
est. The seminal work along these lines is found in 
Refs. 0-0 • This work underlies many of the analy¬ 
ses and exploratory studies of experiments at present 


and future neutrino facilities, including our earlier 

work [THIcj. 

The present paper was undertaken, and used, for 
the purpose of independently confirming the results 
of Refs. 0-H3- We find that the accuracy of the 
expanded oscillation probability is restricted by the 
presence of branch points in the analytic structure of 
the eigenvalues of neutrinos propagating in matter. 
We also show that the regions where the expanded 
results are reliable is different for expansions in a f|J 
and sin 2 $13 SI- We then map out regions where 
the expanded results are reliable by comparing nu¬ 
merical results to the exact results of our Hamilto¬ 
nian formulation. 

Another recent study El takes a complementary 
approach and finds that the predictions of Refs. [7- 
10] can be improved in certain regions using an exact 
evaluation of the integral / Q * rather than the ap¬ 
proximate one found there. It concludes that within 
these regions, predictions of (/i, e) oscillations im¬ 
prove for certain values of the experimental param¬ 
eters. 

The dimensionalities of the neutrino Hamiltonian 
and the parameter space characterizing the mix¬ 
ing of three neutrino pairs are sources of difficulty for 
finding a tractable representation of the oscillation 
probability. The Lagrange interpolation formula Q 
is enormously helpful, providing an exact and for¬ 
mally elegant expression for the exponentiation of 
anixxn matrix. 
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The description of two-flavor neutrino oscillations 
is elementary by comparison. In that case, H v is a 
2 x 2 matrix, and the mixing is described by a single 
real parameter. 

II. NEUTRINO DYNAMICS 

We will be interested in the dynamics of the 
three known neutrinos and their corresponding anti¬ 
neutrinos in matter. This dynamics is determined 
by the time-dependent Schrodinger equation, 

> = H v \v{t) > , (1) 

where the neutrino Hamiltonian, 

H v = H 0v + Hi , (2) 

consists of a piece Hq v describing neutrinos in the 
vacuum and a piece H i describing their interaction 
with matter. 

The solutions of Eq. m may be expressed in terms 
of stationary-state solutions of the eigenvalue (EV) 
equation 

l^mi ^ > •> (A) 

where the label “to” indicates neutrino mass eigen¬ 
states, as distinguished from their flavor states some¬ 
times denoted the label “/”. In operator form, this 
dynamics may be expressed in terms of the time- 
evolution operator S(t',t), which describes com¬ 
pletely the evolution of states from time t to t' and 
also satisfies the time-dependent Schrodinger equa¬ 
tion. 

We will examine neutrinos propagating in a uni¬ 
form medium for interactions constant not only in 
space but also time. Because the Hamiltonian is 
then translationally invariant, attention may be re¬ 
stricted to states, both in the vacuum and in matter, 
characterized by momentum p and therefore having 
the overall r-dependence e lp ' r . In this case expres¬ 
sions may be simplified by suppressing the overall 
plane wave, a convention we adopt. 

For time-independent interactions, S(t',t), 

S(t',t) = e~ lH ^ t '- t) , (4) 

depends on time only through the time difference 
t’ — t. Then, written in terms of the stationary state 
solutions | u m i > of Eq. 0 , 

= ^2 I v mi > e~ lEi (t _t) < v mi \ . ( 5 ) 


With the momentum dependence factored out, 
three basis states | M{i) >, i = (1,2,3) are then re¬ 
quired to describe three neutrinos. The basis states 
correspond to a specific representation, as in descrip¬ 
tions of a spin-1 object. The basis should, of course, 


be orthogonal, 


< M(i)\M(j) > = Sij 

( 6 ) 

and complete, 


J2\M(k) ><M{k )| = 1 . 

(7) 


k 


Once the basis is chosen, wavefunctions for a neu¬ 
trino state are naturally introduced as the compo¬ 
nents of this state in the chosen basis. For example, 
with the eigenstates of Eq. ([3]) expanded in the basis, 

Wmj > = Y^\ M i i ) >m ) > ( 8 ) 

i 

the wave functions of \v m j > would be the set 
to* , i = (1,2,3). With the plane wave factored 
out, the wave function is just a set of three num¬ 
bers. Additionally, introduction of a basis makes it 
possible to represent neutrino states and operators 
such as H„ in matrix form, with each entry in the 
matrix corresponding to a projection of the object 
being described onto the basis. 

In this paper we take the Hamiltonian in Eq. ([2]) to 
be expressed in the standard representation, where 
the mass basis states | M(i) > are taken as the set of 
states that diagonalize the neutrino vacuum Hamil¬ 
tonian H 0v , i.e. | M(i) >= |z^ >= >, 

H 0v \v 0 mi >=E°\v 0 mi > . (9) 

In matrix form 

f E° 0 0 \ 

H 0v = [ 0 El 0 I , (10) 

\ 0 0 E\) 

with the EV’s taken to be ordered E° < E® < E ° as 
in the normal mass hierarchy. In the literature, the 
Hamiltonian is often expressed in a different basis 
obtained by rotating to one in which the complete 
neutrino Hamiltonian is diagonal as in Ref. Q. 

We assume here that that neutrinos and anti¬ 
neutrinos represented by | vf ni > and | r/f vi >, respec¬ 
tively, are the structureless elementary Dirac fields 
of the the Standard Neutrino Model (2]. For this 
reason the theory is invariant under CPT, so the 
mass of an anti-neutrino in the vacuum is the same 
as that for its corresponding neutrino. 
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A. Flavor and Mass States 

Neutrinos are produced and detected in states of 
good flavor, \vji >. The three flavors, electron (e), 
muon (n), and tau (r) correspond, respectively, to 
the index values i = (1,2,3). In the vacuum, each 
flavor state is a specific linear combination of the 
three mass eigenstates \M{i) > of the neutrino vac¬ 
uum Hamiltonian Hq v . This linear combination is 
expressed in terms of the same set of numbers Uij 
for both neutrinos and anti neutrinos 

Wfi > = ^2 u *j\ M {j) > 

3 

\u%> = Y J UiA M {i)> , (ii) 

3 

where Uij are the elements of a unitary operator U, 
the neutrino analog of the familiar CKM matrix. It 
is standard to express Uij in terms of three mixing 
angles ( 0 i 2 , 013 , 023 ) and a phase S cp characterizing 
CP violation, 

/ C 12 C 13 S 12 C 13 s 13 e~ lS ^ \ 

I U 21 U 22 S 23 C 13 I , ( 12 ) 

\ U 3 \ U32 C23C13 J 

where 

U 21 = —S 12 C 23 — ci2S23Si3e i<5op 

U 22 = C12C23 — Sl 2 S 23 Sl 3 e* l5 ' =p 

U31 = S12S23 — Ci2C2 3 Si 3 e l5cv 

U32 = — C12S23 — Si 2 C 23 Sl 3 e l,5oP • ( 13 ) 

We use here the standard abbreviation S 12 = sin0i 2 , 
C 12 = cos 0i 2 , etc. The parameters 0 and S cp are 
determined from experiment. 

Because Uij —> U*j with S cp —> —S cp it follows 
that the relationship in Eq. m between flavor and 
mass states for anti-neutrinos and neutrinos in the 
vacuum is equivalent to S cp 0 —S cp . 


B. Neutrino Interacting Hamiltonian 

The interaction Hi, determined by taking the elec¬ 
tron flavor states scattering from the electrons of the 
medium to mediate the interaction, is then expressed 
as an operator in the standard representation, 



with V = P\f2Gpn e and n e the electron number 
density in matter. 

For electrically neutral matter consisting of pro¬ 
tons, neutrons, and electrons, the electron density 
n e is the same as the proton density n p , 

Tl e = Tip 

= RN , (15) 

where N = n n +n p is the average total nucleon num¬ 
ber density and R = n p /N is the average proton- 
nucleon ratio. In the earth’s mantle, the domi¬ 
nant constituents of matter are the light elements 
so R « 1/2; in the surface of a neutron star R << 1. 
Matrix elements of Hi are thus 

< M(k)\Hi\M(k') > = U{ k VUi k , . (16) 

Matrix elements of Hi are thus 

< M{k)\Hi\M{k') > = U*i k VUi k , . (17) 

C. Dimensionless variables 

The results are most naturally expressed in di¬ 
mensionless variables. We first take advantage of the 
global phase invariance to express all energies rela¬ 
tive to the vacuum EV E® of the same momentum. 
We indicate that a quantity is expressed relative to 
E® by placing a bar over it, e.g., 

= E° - E° . (18) 


We follow the same convention for the Hamiltonian, 


H v = H V - 1 E° , 

(19) 

so the EV equation Eq. ([3]) becomes 


(Hq v + i7i)|zz m j > = E^v^i >, 

( 20 ) 

where 


Hov = H 0v — IE® . 

( 21 ) 


Then, to express the theory in dimensionless vari¬ 
ables we divide all energies, including the Hamil¬ 
tonian, by E 3 = £3 — E®. The stationary-states 
\v m i > are also be determined from the dimension¬ 
less Hamiltonian H„, 

Hu = Hov + Hi , ( 22 ) 

i.e., from the solutions of 

Hu\Umi ^ = Ei^Umi , 


H 1 = U 


V 0 0 
0 0 0 
0 0 0 


(14) 


(23) 
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where the “hat” placed over a quantity indicates it 
is dimensionless. Thus 


A o. = §jf, (24) 

and Hi is obtained from Hi by replacing 

(25) 

The quantity A is the same as that defined in 
Refs. I4i. I7H11|. The connection of the Hamiltonian 
H„ to the full Hamiltonian H v = Hq v + Hi is then 

H v = 1 E°i + E$H„ . (26) 


D. Neutrino Vacuum Hamiltonian Hq v 


The case of main interest for many situations is 
the ultra-relativistic limit, \p\ » m 2 (we take the 
speed of light c = 1). For ultra-relativistic neutrinos 
in the laboratory frame, the energy of a neutrino in 
the vacuum becomes 

9 

TT) * 

« bl + ^ , (27) 

where rrii is its mass the vacuum. Similarly, Ei ap¬ 
pearing in Eq. 0 may be written 

M 2 

Ei « bl + > (28) 


where Mi is its mass in the medium. Thus, in this 
limit, 


and 


with 


Ei —> 


M 2 — to 2 
2 2~ 

m 3 ~ m l 


/ 0 0 0 \ 
Hoy —^ [ 0 a 0 j 

V° 0 !/ 


2 2 

777,2 — 777, J 

2 2 ' 

77^3 — 777,J 


(29) 


(30) 


(31) 


In this limit, the distance L from the source to the 
detector corresponding to S(t',t ) in Eq. 0 is 


(32) 


The time-evolution operator, Eq. 0 , expressed in 
dimensionless variables is, 


S(L) = 

_ e 2iE 1 AL e —2iH^Ar, 


(33) 


where H u is given in Eq. (1261) . and where 


L(m 2 - ml) 
L ~ 4 E 


(34) 


[The similar quantity as defined in Ref. 7| is 
exactly one-half of that appearing in Eq. (1M1) .] 


E. Neutrino Mass Eigenvalues 

The neutrino mass eigenstates in a medium are 
solutions to the EV equation for H u , Eq. ([10. In 
many familiar formulations the full solutions, 
including both the eigenstates > and EV’s Ei, 
are required to find the neutrino oscillation proba¬ 
bilities. 


1. Diagonalization of Neutrino Hamiltonian 


The energies Ei are solutions of the cubic equa¬ 
tion jl 2 j 

Ei + aEi + bEi + c = 0 , (35) 

where 

cl = —(1 —F cr —f— H) 

b = a + A cos 2 0 i 3 + AaC^ 

C = —Aol cos 2 012 cos 2 013 . (36) 

We have expressed b in terms of a frequently recur¬ 
ring combination of mixing angles, 

= cos 2 0 i 2 ± sin 2 0 i 2 sin 2 0 i3 . (37) 

Note that the mass eigenstate energies are indepen¬ 
dent of S cp and 023 for both neutrinos and antineu¬ 
trinos. 

The solutions of Eq. (1551) are expressed conve¬ 
niently in terms of the quantity d, 

d = tp + y /1 p 2 — 4y 3 
7 = a 2 — 36 
ip = a 3 — 27c — 3aj . 


L = t' -t . 


(38) 
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These solutions are real when 

|rfi/3|2 = 2 2 / 3 7 > 0 , (39) 

which requires 

V' 2 < 4y 3 , (40) 

and, thus, that d be complex. Because having real 
energies is required by Hernriticity of the neutrino 
Hamiltonian, Eqs. (1391401) amount to conditions on 
all parameter sets in terms of which H „ is defined. 
We find 

Ei = ~ - g. 21/3 (V^ Im[d 1/3 } + Re[d 1/3 ]) 

Ei = -f + j^j-^VZImid 1 ' 3 } - Reid 1 ' 3 }) 
a 2 2//3 

E 3 = -- + — Re[d 1/3 } . (41) 

The masses are ordered so that m 3 > m 2 > m 1 . 
Because EV do not cross, E 3 > E 2 > E\ for all 
|A|. A simple constraint among Ei is found from 
the trace of Eq. (1^1) , 

TrH v = E\ + E 2 + E 3 
= TrH 0 v + TrHi 
= 1 + a + A = —a . (42) 

2. Using neutrino mass eigenvalues in our 
Hamiltonian Formulation 

In our formulation, the entire dependence of the 
time evolution operator on the neutrino eigenvalues 
occurs through three eigenvalue combinations, 

A Eiii = Ep — Ep> 

E«< = Ei + Epi 

n«/ = EpEpi , (43) 

with £ > l' (and powers thereof). We denote such 
quantities using a bracket notation, For example, 

AE[1] = E 3 -E 2 
AE[2] = E 3 -E 1 

AE[3] = E 2 -E 1 , (44) 

in the case of A E. We will generally use this bracket 
notation also for other quantities in our formulation 
that depend on two indices (£,£'), such as and 

n 


An expression for E[^], 

E [£] = -a-E e , (45) 

follows from Eq. (1421) . An equivalent expression for 
n[£] in terms of Ep is found by subtracting Eq. (1351) 
for Ep and that for Ept and dividing through by 
A Eppt. We find 

0 = (Ep + EpEpf + Epi ) 

+ a(Ep + Epi ) + b 

= (Ep + Epi) 2 - EpEpi 

+ a(Ep + Epi) + b , (46) 

giving 

T,[£} 2 - n[t\ + aY,[e] + b = 0 . (47) 

Then, using Eq. (Il5l) . 

n[£] = Tip(Yip + a) + b 

= b + aEp + Ep . (48) 

Finally, having observed that powers of the quan¬ 
tities given in Eq. IS will appear in various expres¬ 
sions, we note that n[£] p and E [£] 9 with p > 2 and 

~ n 

q > 3 involve linear combinations of eigenvalues Ep 
with powers n > 3. Such terms are equivalently rep¬ 
resented by a linear combination of three terms, one 

- 2 ~ 

proportional to Ep , one proportional to Ep , and one 
independent of Ep, obtained by using the equation 
of motion repeatedly. We later make use of this fact 
to simplify various expressions. 

III. THE S-MATRIX IN OUR 
HAMILTONIAN FORMULATION 

The probability V(u a —► up,) for neutrinos to oscil¬ 
late from the initial state of flavor a to a final state 
of flavor b is found from the time-evolution operator 
S(t',t) as 

V(v a ^v b ) = \S ab (t',t)\ 2 

= P ab (t' - t) , (49) 

where 

\S a \t',t)\ 2 = \<v%\S(t',t)\v%>\ 2 . (50) 

We accordingly determine here P ab (t' — t) from 
S(t',t) defined in Eq. £H). 
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In this section we review the formulation of neu¬ 
trino oscillations based on the Lagrange interpo¬ 
lation formula as used in Ref. [3j. This formula¬ 
tion leads to exact, closed-form expressions for the 
time-evolution operator and the partial oscillation 
probabilities that are linear combinations of spher¬ 
ical Bessel functions in the eigenvalue differences 
whose coefficients are polynomials in the neutrino 
CKM matrix elements, the neutrino mass differences 
squared, the strength of the neutrino interaction 
with matter, and the neutrino mass eigenvalues in 
matter. We are led quite naturally to such expres¬ 
sions for all the partial oscillation probabilities in 
terms of these basic quantities. The numerical re¬ 
sults given later in this paper are based on this for¬ 
mulation. 


A. Time-Evolution Operator 


The overall phase in Eq. (1531) does not contribute 
to \S ab (t', t)\ 2 , so for the purpose of calculating the 
oscillation probability, we may take 

S(L) -> e ~ ik ^ L . (51) 

Then, with neutrinos created and detected in flavor 
states, which are coherent linear combinations of the 
neutrino vacuum mass eigenstates given in Eq. ©, 

Wfa> = Y, U *ai\ M (ti> , ( 52 ) 

j 

we see that the mass eigenstate components of 
the flavor states contribute coherently to the time- 
evolution operator. Thus, 

/ ,,0 \ -iH„A L I 0 v 

< i/f b \e \Vf a > 

= <M(b)\Ue- iS “ AL Ui\M{a) > . (53) 

This coherence leads to the oscillation phenomenon. 

The elegant formulae for S(L) = are 

obtained from the Lagrange interpolation formula, 
Eqs. (9,11) of Ref. Q, 

Ue -iH»A LU -1 = ^2 F(exp -iE e A L , ( 54 ) 

i 


where T = L = t' — t and 


Ft = 


UHyU- 1 - lEj 

Et — Ej 


(55) 


For three neutrinos, the sum in Eq. (1541) runs over 
three values of £ and the product in Eq. (l55l) over 
two values of j. 


Using the convention that O ab , written without 
parentheses enclosing ab , denotes the matrix ele¬ 
ments of the operator O, 

O ab = < M(b)\0\M(a) > , (56) 

the matrix elements F^ b of Ft given in Eq. (l55l) may 
be compactly written 


rpcib 


where 0, 


< M(b)\W[£]\M(a) > 

m 


(57) 


W[ 1] = {UH U U~ 1 - 1 E^UHvU- 1 - 1 E 2 ) 
W[2\ = (UH„U~ l - 1 E 3 )(UH v rf - lEi) 

W[3] = {UHvU- 1 - 1 E 2 ) 

x (UH„U ] - lEx) (58) 

and 

D[1] = (E 3 -E 1 )(E 2 -E 1 ) 

D[ 2] = (Ei - E 2 )(E 3 — E 2 ) 

D[3] = (E 1 -E 3 )(E 2 -E 3 ) . (59) 

Equations (I58I59D use the same bracket notation in¬ 
troduced in Eq. (1551) . The result in Eqs. (I53I54I55I) 
is immediately verified to be correct by inserting a 
complete set of intermediate eigenstates of H v in 
Eq. (1551) . 

It follows from the unitarity of U that W[£] is 
Hermitian, 


W [£] f = W[i] (60) 

and that the two factors in Eqs. (1551) commute with 
each other. We find from Eq. (l58l) that 


TrW[£} = D[£] . (61) 

Equation m establishes the reflection symmetry 
Ff* = Fj! a . (62) 


Explicit expressions for W[£] are easily found in 

terms of H v . The entire dependence of W[£] on S cp 
occurs through three operators independent of 5 cp , 
W 0 [£\, W cos [£} and W sin [£], 


W[£\ = W 0 [£] + cos6 cp W cos [£] 

+ i sin 5 cp W sin [£] , (63) 


with Wo[^], W cos [£] and W S i n [£] real and indepen¬ 
dent of S cp . Details are given in Appendix 1X1 
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B. Total Oscillation Probability 

Expressions for V{v a —> v b ) may be obtained di¬ 
rectly from S{L), 

V(v a ^v h ) = \ S ab (t',t)\ 2 

= Re[S ab (L)] 2 + Im[S ab (L)] 2 . (64) 

Convenient expressions for Re[S ab (L)] and 
Im[S ab (L)] defined by Eq. (15-0) are presented 
in Appendix [A] In our Hamiltonian formulation, the 
dependence of S(T) on the CP violating phase S cp 
is very simple and follows from Eqs. (1541631) noting 

a ab ±. 

that F“ b = W [£]/D[£], Eq. (1571) . We thus find, 


^ ab 

Re[S ab (t\ t)} = 6 ab - 2 V sin 2 E t A L 

V m 

± ab 

— 2 COS S cp OS [l] s j n 2 

i m 

± ab 

_ W • \£] ± 

sin 5 cp f"l sin 2E t - A l (>65) 


m 


and 


~ ab 

Im\S ab (t ', t)} = -2 sin S cp V W sin 2 E t A L 

V m 

n ab 

_y^^oM sin2 E eAL 
t m 

- ab 

E W \£] - 

— f os sin2F^A^ (66) 

t m 

where is defined in Eq. (l34l) . Approximate ex¬ 
pressions for P(v a —> Vb) in terms of the parameters 
of H„ were obtained from S(L) in Refs. [H, @| by an 
expansion in sin(?i 3 - 


C. Partial Oscillation Probabilities 

Using somewhat different techniques, the oscilla¬ 
tion probability may be expressed through a set of 
functions that express how V(y a —> v b ) = P ab de¬ 
pends on the CP violating phase S cp In our 
Hamiltonian formulation there are four such terms, 

P ab = S(a, b ) + P 0 “ b + P s t s + P c ab s5 


with P s “ b 5 linear in sin<5 cp , P c a 0 b S(5 linear in cos<5 cp , 
P^os 2 s quadratic in cosS cp , and Pq 6 independent of 
S cp ■ Although only the overall oscillation probabil¬ 
ity is a true probability, guaranteed to be strictly 
positive everywhere, we find it convenient to refer 
to these four terms as “partial oscillation probabili¬ 
ties.” Approximate expressions for the partial oscil¬ 
lation probabilities expanded in the small parameter 
a of the SNM in Ref. [||. 

Obtaining expressions for the partial oscillation 
probabilities from the time-evolution operator re¬ 
quires additional analysis, given in Appendix [B] In 
terms of spherical Bessel functions, we find there, 


PsinS (Al i A) = sin 5 C p^ £(- 1)V 


D 

x io(2A[f]) , 


ab 

sin 




( 68 ) 


where w“ b n [£], and therefore Pg b n g , is anti-symmetric 
under a b. The other three partial oscillation 
probabilities are individually symmetric under a 
b. We find 

4A 2 

P coss( a l,A) = - cos $ cp — 

D V 

x AE[£}j 2 (A[£}) 

A A 2 

Pcos 2 s(A L , A) = -cos 2 4p^E(-l)'<« 

D p 


b 

cos 2 


x AE[£]j 2 (A[£\) 
4 A 2 

P 0 “ b (A i ,A) = -^^(-l) 

D V 

x AE[i\j 2 {A[t\) , 


\£ n .,ab r 


W , 


0 


(69) 


P 


ab 


S 2 8 5 


where all sums run over £ = 1 ,2 ,3, A [t\ is defined 
as 

A[£] = AE[£\A L , (70) 

with A E[£] defined in Eq. (l44l) . D is defined as 

D = AF[1]AF[2]AF[3] , (71) 

and the matrix elements w“ b \£] are given in terms of 
the mixing angles and A in Appendix [B] Since we 
order the energies so that E 3 > E 2 > E 3: AE[£] as 
well as D are all positive. 

We begin our derivation with the expression for 
the oscillation probability written in terms of S(T ), 
Eq. 454}, 

V(y a v h ) = P ab 


(67) 














= | < M(b)\Ue- iH ^'- t] U ] \M{a) > | 2 

= J2 F M' ex P~ i{Be ~ Be ' )L ■ 

w 


Here F° b , is defined as 


rpcib _ rpab rpab* 

*1’ 


D[£]b[£') 


with D[£] given in Eq. (1591) . and 

10 $ = < M(b)\W[£]\M(a) > 
x < M(b)\W[f]\M(a) >* . 


(72) 


(73) 


(74) 


All results needed for determining exact, closed-form 
expressions for the partial oscillation probabilities 
are found in Appendix E 

As we explain in Appendix E w $, may be ex¬ 
pressed through four operators, 

w w = w o,it< + cos S c P w“os w + c °s 2 S cp 

X KoaHC +*sin 5 C pwf inU , . (75) 


found from the decomposition of W[£] given in 
Eq. (l63l) . The quantities w° ,b [£] appearing in 
Eqs. (|68I69I) are the same as wf b e , written in the 
streamlined notation, 


<'[!]=<«' for (44) = (3,2) 
wf[2]= w% for (*,0 = (3,1) 
wf[3]= w% for (44) = (2,1), (76) 

which takes advantage o£ £ > £'. 

An analytic expression for Wg b n [£], 


<1 M = K sin 29 23 AE[£]e a s b D , (77) 

follows directly from Eq. (1951) . Here, e s i n is the anti¬ 
symmetric matrix 


0 1 —1 \ 

-10 1 (78) 

1 - 10 / 


and 


K = - 


a(l — a) 


x cos 0 i 3 sin 2012 sin 2013 . 


(79) 


Equation fzzD is one of the more striking results. 
Analytic formulae for the other w°’ b [£] follow from 
Eq. (1571) . These are all given in Appendix E where 
they are expressed in terms of the parameters spec¬ 
ifying Hq v and H\ = U~ 1 VU. These are the same 
parameters defining SNM. 


It follows from Eq. m that the coefficients of 
AE[£]j 0 (2A[£}) in Eq. (l 68 l) are all proportional, lead¬ 
ing to a simple expression for Pg b n g, 

P sLs(Al, A) = A|a(l - a) sin<5 cp e“ b n 

x cos 0 i 3 sin 2012 sin 2013 sin 20 2 3 
x jo (A [3] )jo (A [2]) jo (A [1]) . (80) 

Note that P^ s is antisymmetric under a -O- b. 

Analytic expressions for all other partial oscilla¬ 
tion probabilities follow from Eq. C 2 } using Eq. Cal), 
which expresses F^ b , in terms of In this fash¬ 

ion, these partial oscillation probabilities are also ex¬ 
pressed in terms of the parameters of the SNM and 
the neutrino eigenvalues, 

The usefulness of the partial oscillation probabil¬ 
ities can be seen as follows. It is a general result 
that the exchange of initial and final states in the 
oscillation probability or neutrinos (antineutrinos) 
is equivalent to letting S cp —> —5 cp . Thus, the re¬ 
sult for the inverse reaction V(v b —> v a ) is found by 
exchanging (a, b) in Eq. (1671) . Since P£ b g is antisym¬ 
metric under the exchange of (a, b), and P§ b , P c “ b S(5 
and P ab 2 g symmetric, it follows that P ba is given 
by 

V[y b u a ) = 8(a, b ) + P 0 “ b - Pf nS + P c ab sS 

+ P cos 2 <5 ■ (81) 

In analogy to Eq. (EH), we may express the oscil¬ 
lation probability for antineutrinos as 

V(y a -»■ D b ) = P ab 

= 6(a, b) + pQ b + P s “n s + PcosS + P cos 2 5 > (82) 

where the bared probabilities for anti neutrinos are 
obtained from the unbarred for neutrinos by replac¬ 
ing S cp —► —5 cp and A —► —A. Because the energies 
of antineutrinos are different from those of the neu¬ 
trinos in matter, we can expect P ab p P ab in this 
situation. 

Again applying the rule that exchange of initial 
and final states is accomplished by letting S cp —> 
—5 cp , the oscillation probability V(v b -* D a ) is ex¬ 
pressed in terms of the same four quantities, 

V(pb -t v a ) 

= 6(a, b) + P 0 “ b - P s f aS + P c ts + P c t* 5 ■ (83) 

It is worth noting that the entire dependence of 
the oscillation probabilities given in Eqs. (I68I69D on 
the neutrino beam energy E, the baseline A, and the 
medium properties occurs through the variables A l 
and A defined in Eqs. (1341251) . respectively. Since 
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we will be most interested in how the neutrino os¬ 
cillation probability depends on the beam energy, 
baseline, and medium properties, the partial oscilla¬ 
tion probabilities have been expressed as functions 
of A l and A. 

Because the vacuum result is also easily obtained 
by less sophisticated arguments, the vacuum limit 
provides an opportunity to verify our Hamiltonian 
formulation in a well-known special case. 


D. Special Cases 

We next examine special cases chosen to highlight 
specific aspects of our formulation. Because these 
cases are well known in other contexts, they pro¬ 
vide useful cross checks. We first consider two-flavor 
mixing and then the vacuum limit. 


1. Two Flavor Mixing 


with Al given in Eq. (1M1) . The two-flavor oscillation 
probability P 2 12 (Al, A) is then 


p$(A L ,A) = \st}f 


sin 2 20 


sin 2 A L cj) , 


(90) 


which we will next compare to our Hamiltonian for¬ 
mulation. To find S), 2 =< M(2)\S[L]\M(1) >, it is 
easy to overlook that the 2x2 matrix H v in Eq. (1M1) 
is not diagonal in the flavor basis . 

The oscillation probability in the three neutrino 
mixing of our Hamiltonian formulation is matched 
to the two-flavor case just discussed by setting two 
of the mixing angles in Eq. (fT^l) to zero, say 9 1 3 —>• 0 
and 023 — > 0 and identifying 9 with 9 12 . Accordingly, 
we find, 


/ cos 9 sin 9 0 \ 

U = [ — sin0 cos 9 0 j . (91) 

V 0 01 ) 


Expressions for two-neutrino oscillations are easily 
obtained from the Hamiltonian, 

H v = H 0 u + U-'H^U , (84) 

where 

( ° J ) , (85) 

A-(oS)- < 86 > 

The standard mixing matrix U is 

/ c osfl S i„0\ 

— sin 9 cos 9 J v ' 

The neutrino Hamiltonian is easily diagonalized to 
find the medium-modified two-flavor neutrino eigen¬ 
values, 

E 2 = \ (l + A + 

i 1 = i(l + i-<A) , (88) 

where </> = y/1 + A 2 — 2 A cos 29. 

Expressed in dimensionless variables, the time- 
evolution operator S ab for the transition a —> b be¬ 
comes 

S ab = < M(b)\Uexp- 2lR ^ AL U~ 1 \M(a) > ,(89) 


This is easily recognized as a two-flavor mixing ma¬ 
trix by noting that one of the three neutrinos does 
not mix with the other two. We see that it not only 
depends on just one mixing angle 0 , but also that the 
dependence on the CP violating phase <5 has dropped 
out. 

Natural choices for a 3 x 3 two-flavor neutrino vac¬ 
uum Hamiltonian and interaction are, 

/ 0 0 0 \ 

4= 010 , (92) 

\ 0 0 0 J 

and 

^ (A 0 0 \ 

H 1 = 0 0 0 . (93) 

\o 0 0 / 

When the 3x3 neutrino Hamiltonian H v = Hqv + 
U~ 1 HivU is diagonalized, two of the three medium- 
modified neutrino eigenvalues, 

4 = \ (1 + A + </>) 

4 = ^ (i + A - 0 ) 

4=0. (94) 

are identical to those found above in the two-flavor 
case. Note that the eigenvalues have no branch 
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points for A on the real axis. From these eigenvalues, 
we find the differences, 

AE[3\ = Ei- Ex 

4( 1 + M 

AE[2] =E 3 -Ex 

= l( i + i + 0) 

AP[1] = E a - E 2 = (j) . (95) 

Finally, consider the oscillation probability in our 
Hamiltonian formulation. Because P^g, P c “ b s< 5 , and 
Pcos 2 s are proportional to sin 29 23 (and 9 23 has been 
set to 0), these terms do not contribute for two-flavor 
mixing. Thus, 

P ab (A L ,A)=6(a,b) + PS\A L ,A) (96) 

with Pq 6 (Al, A) taken from Eq. (1691) . We find after 
a straightforward calculation, 

wl 1 [e] = wl 2 [£] 

wl 2 [e] = -wl 2 [£\, (97) 

where, 

Wq 2 [^] = {A cos 2 8 — (1 + A)Ei 

+ E t ) . (98) 

Taking Eg from Eq. (PI) . we evaluate 

A cos 2 9 - (1 + A)E e + E e , (99) 

to find 

®o 2 [1 ]=«'o 2 [2 ]=0 , (100) 

and 

w q 2 P] = A cos 2 dw^ 

A 

= — cos 2 9 sin 2 29 . (101) 

4 

I 


Thus, for all two-flavor transitions, 

^ 4 A 2 

P ab (A l , A) = 8(a, b) + -^< b [3]AP[3] 

x j 2 0 (AE[3]A l ) , (102) 

where D is from Eq. (11031) . 

D = AP[1]AP[2]AP[3] 

= ^|^((l + i) 2 -l-i 2 + 2 icos 26 i) 

= icos 2 6 »AP[l] . (103) 

To compare our result to to the well-known two- 
flavor oscillation probability given in Eq. (l90l) . we 
evaluate Eq. (11021) for 1 — > 2 transitions obtaining, 

P 12 (A l , A) = sin 2 A L (j) . (104) 

<r 

As expected, this is in complete agreement with 
Eq. m . 

2. Vacuum Oscillation Probability 

In our Hamiltonian formulation, an expression 
for the time-evolution operator in the vacuum limit 
A — > 0 is found from Eqs. (ESI), (ED, (53, and (EH 

< M(b)\S°(t',t)\M(a) >= S 0ab {t' - t ) 

= E <M(5)|W°[£]|M( Q )> ex ^ 

t D°[£] 

Here, E® and D° [£} are the vacuum values of Eg and 

D\£\, respectively. In the vacuum, of course, there 
is no distinction between the energy differences for 
neutrinos and anti-neutrinos. 

Evaluating Eq. (11051) by inserting a complete set 

a 0 

of states intermediate states |n >< n\ inside W [l], 
we arrive at 

(106) 


< M(b)\S°{t',t)\M(a) > = EE U bl 


-0 -0 -0 ~0 
(E n -E a )(E n -E b ) 

D°[£] 


U* e 

u na c 


-iEg {t' ~t) 


The indices (a, b, n ) run over all permutations of the 


three integers (1,2,3), from which it follows that 
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-0 ~0 -0 -0 _~ 

( E n - E a[e])( E n ~ E b[i\) = D °l 4 and, consequently, 

< M(b)\S°(t',t)\M(a) > 

= Y,U bn e- iE W-^U* na . (107) 

n 

This well-known vacuum limit also follows di¬ 
rectly from more elementary considerations, using 

Eqs. mm, 

< Vf b \ e ~ iH ° v ^ _t ) |i/° a > 


= Y, U ^~ iE ° n(t '~ t)U L ■ ( 108 ) 


The equality of Eqs. (110711081) verifies our exact 
Hamiltonian formulation in the vacuum limit. 

Finally, we find that Ffy appearing Eq. m takes 
a familiar form in the vacuum. Using Eq. (11071) . or 

Eq. cna, 




) = (E Uwe-^W-Vu^ (^E Ubfye-mdt'-Dup^ 


Expressing this in terms of Freund’s 

and using Eq. (EZD, 


It' Hi’ 

Ju- = u; e u ai u: e ,u w 


TT'O rpO 

J-J 0 J-Jpr — 


2 E 


we find that Eq. (11091) may be written, 

V(v° a -+ v° b ) = 5 ab - 4 E ReJtj sin 2 ((£° - E° e ,)L ) - 2 E sin(2(U° - E°)L) . 


(109) 

( 110 ) 

( 111 ) 

( 112 ) 


i>j 


i>j 


Comparing Eqs. (11121) and d72l) . it is clear that 
F°j, —> jfy* in the vacuum limit. 


IV. THE STANDARD NEUTRINO MODEL 

We adopt the Standard Neutrino Model (2J as 
our description of neutrino physics. The parame¬ 
ters defining the model include a (dimensionless) in¬ 
teraction strength A of neutrinos and anti-neutrinos 
with matter, the three neutrino mass differences, the 
three mixing angles, and the CP -violating phase. 
Most of the parameters of the SNM are consistent 
with global fits to neutrino oscillation data with rel¬ 
atively good precision. 

The neutrino mass differences of the SNM are 
taken to be [zHU 

2 2 _ r 2 

m 2 — m 1 = om 2 i 

= 7.6 x 10 " 5 eV 2 

_2 _2 _ e_2 

ra 3 — m 1 = om 31 


= 2.4 x 10 " 3 eV 2 , (113) 

corresponding to 

= i 

a ~ Sm 31 

= 3.17 x 10~ 2 . (114) 

In Ref. @-S3|, <5 = (m| — m\)/(2E) and A = (m§ — 
mf)/(2E) whereas in Ref. [5J|, A = {m%—m\)L/(AE) 
and A = m\ — m\. 

The mixing angle 023, 

023 = tt/4 , (115) 

is the best-fit value from Ref. j43j], and 0 i 2 , 

0 12 = tt/5.4 , (116) 

is consistent with the recent analysis of Ref. 0 . 
The mixing angle 0 i 3 is known to be small ( 0 i 3 < 
0.18 at the 95% confidence level) but its precise value 
is uncertain. A very recent result from the Daya 
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Bay project M is sin 0 i 3 « 0.15, which we adopt to 
determine our value for $ 13 , 

0 13 = 0.151. (117) 

This fixes R p = sin 2 9\ 3 /a ft 0.711. The CP violat¬ 
ing phase 5 cp is not known at all and will one of the 
major interests at future neutrino facilities. 

Parametrizing the interaction strength A, 
Eq. E51) . we find 

A = ±6.50 10” 2 i? E [GeV] p [gm/cm 3 ] , (118) 

with E[GeV] being the neutrino beam energy E (in 
GeV) and p[gm/cm 3 ], is the average total density 
(in gm/cm 3 ) of matter through which the neutrino 
beam passes on its way to the detector (the matter 
having average proton-nucleon ratio R). For our cal¬ 
culations we are interested in experiments close to 
the earth’s surface, so we take 

p[gm/cm 3 ] = p 0 

= 3 , (119) 

the approximate mean density of the earth’s mantle. 

In the SNM, A^, defined in Eq. (i34l) . may be 
parametrized in the high-energy limit as 

AiasM 5 xl 0 "l^j' (120 > 

Here L[Km] is the baseline and A [GeV] is the neu¬ 
trino beam energy. 

V. EIGENVALUE EXPANSIONS IN THE 
SNM 

The fact that a and sin $13 are naturally small 
in the SNM commonly motivates approximation 
schemes [11, 0 based on first-order Taylor se¬ 
ries expansions in one of these small parameters, 
(where ^stands for a or sin 2 $ 13 ). For example, in 
Refs. [ 5 !, ^ the oscillation probability is expanded in 
sin 2 0i3. Reference [ 4 ] makes use of an expansion in 
the small parameter a. 

Although these expansions may be used effectively 
to simplify the theory, they come at a price [H. This 
price is that neither expansion gives accurate repre¬ 
sentations for all values of the interaction strength 
A, including values in some regions of critical impor¬ 
tance. 

In our subsequent work [0, we will consider sim¬ 
plifying the oscillation probability using the same ex¬ 
pansions, but they will be used somewhat differently, 


in two stages. In the first stage, the eigenvalues are 
expanded, as below. In the second, the expansions 
will be used to simplify the oscillation probability. In 
this work, the eigenvalue expansion will, of course be 
made before introducing R p . For the expansion of 
the oscillation probabilities, R p may be introduced 
before the expansion is made because of its simpler 
analytic structure. 


A. Analytic Structure of Eigenvalues 

The key for identifying which of the ^'-expansions 
might be appropriate over specifics ranges of A is 
revealed by the analytic structure of the eigenvalue 
Ee(^). It is particularly important to identify the 
locations of its branch points when £' = 0. Branch 

points identify where a series expansion of EAC), 
or a function of it such as < F(H„) >, would not 
converge. 

We have seen that £/(£) depends on £ entirely 
through the two functions (U(£), V(£)), 

Et(Z) = Et(U(£),V{£)) , (121) 

where 

Z7(0 = Re[d 1/3 } 

V(0 = Im[d 1/3 ] , (122) 

with d defined in Eq. (l38l) . 

Recalling that 4y 3 — ip 2 > 0, we conveniently write 

U {0 = ±*\/ 47 3 - i > 2 ) 1/3 

+ (i/j - i\J 4y 3 - ■i/' 2 ) 1 / 3 , (123) 

and 

^(0 = + iV^T) 1/3 

— ('!/ — *\/4y 3 — 1 Z 2 ) 1 / 3 . (124) 

Branch points clearly occur for parameter values 
satisfying 

V’ 2 (C B )=4 7 3 (e B ) . (125) 

The locations of the branch points of EAA) are 
found as follows for the two small-parameter choices 
& 

When £' = a, branch points as a function of the 
other parameters of the Standard Neutrino Model 
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are found from Eq. (11251) with 'i/>(£ s ) and 7(£ S ) eval¬ 
uated in terms of 

a(^) = -(l+i) 
b(£ B ) = Acos 2 0i 3 

c(C B ) = 0 , (126) 

where t; B = £| a =o- The only real solution of 
Eq. (11251) is 


A = 0 . (127) 

When £' = sin 2 0i 3 , branch points as a function 
of the other parameters of the SNM are found from 
Eq. (11251) with '0(£ s ) and q(£ s ) evaluated in terms 
of 


a(£ s ) = -(1 + A + a) 

5(£ s ) = a + Aa cos 0i2 + A 

c(£ s ) = — Aa cos 2 0 i 2 , (128) 

where t; B = £|e 13= o. The only real solution of 
Eq. (11251) is 


1 — a cos 2 0 i 2 

« 1 . (129) 


yaRp 
2 C B 


(2 — 4A 0 Aq — 2 ol T 2Aqol) 


— — (A + a + Cg) 


OL R>p A 'J . 

(2 — Aq — a + Cg) 


2(1 - y)C e 
yaR p 


2 Cg 


(2 — 4Aq -j- Aq — 2 oi 2Aoa) 


ve _ i i aRpA- 

3 ~ + t^7’ 


(131) 


where y = A/Aq, with Aq = A e 0 the location of the 
branch point for the sin 2 6\q expansion and 


Cg = \JA 2 + a 2 - 2 ia cos 20 i 2 . (132) 

Note that E® = Ef Above the branch 

point, y > 1, and E% exchange roles. 

The lack of convergence at the branch point is 
manifest here by inspection, he., by the appearance 
of simple pole at y = 1. Although Eq. (11311) would 
suggest a pole in all three eigenvalues, E ® is rather 
accurate for all values of y. The absence of a singu¬ 
larity in Ei can be confirmed by a simple calculation 
that shows the coefficient of (1 — j/) -1 vanishes at 

y = i- 


B. The First-order Taylor Expansions 

Defining E e to be the first two terms of the Taylor 
series for Eg expanded about £ = 0, we find 

E e = En| 5=0 +C-^7-|{=o j (130) 

where E( and dE^/d^ are easily obtained from 
Eqs. (1351361) and their derivatives. The results for 
£ = sin 2 0 i 3 and £ = a found in this way are given 
immediately below. 

1 . f = sin 2 0i3 

Applying Eq. (11301) we find for £ = sin 2 0i 3 and 
below the corresponding branch point, 


2. £ = a 

Taking R p = sin 2 0i 3 /a, we may write cos 2 0i 3 = 
1 — aR p and cos20i 3 = 1 — 2 aR p . Then, for £ = a 
and y > 0 (above its corresponding branch point), 

Ei = a cos 0 2 2 

E$ = l{l+A-d a ) 

sin0? 2 1 — A( 1 — 2aR p ). 

+ a 0 (1 +-^-—) 

2 

Eq = -(1 + A + C a ) 

+ _ 1 -ia- 2 ^)) _ (133) 

2 C a 

where 


E{ = -(A + a- Cg) 


+ 


aRp A 
2(1 - y)C B 


(2 — Aq — a — Cg) 


C a = y/( 1 - A) 2 + 4 olR p A . (134) 

Note that Eg = Eg • These expressions are 
identical to Eqs. (18,19) of Freund [||. 
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In contrast to £ = sin 2 0i 3 , the lack of convergence 
of the a expansion is not obvious from a casual ex¬ 
amination of Eq. (I133D . Numerical comparison to 

the exact result confirms that £7“ and E% are poor 
representations of the corresponding exact results in 
the vicinity of the branch point at y = 0. However, 
no evidence of the branch point is apparent in , 
which is rather accurate for all values of y. Across 
the branch point at y = 0, E“ and E% exchange 
roles. 


C. Branch Points and Resonances 

Resonances are heralded by the appearance of 
minima in the EV differences. Two well-known neu¬ 
trino resonances occur in the SNM. One of these, 
the solar resonance, is found for relatively weak in¬ 
teraction strengths A ss a. The other, the atmo¬ 
spheric resonance, occurs for stronger interactions, 
A sa cos 2013 . The solar resonance occurs very close 
to the branch point identified with the a expan¬ 
sion, and the atmospheric resonance very close to 
the branch point identified with the sin 2 0i 3 expan¬ 
sion. 

Using Eq. (111 8 [) . the solar resonance appears for 
neutrinos of energy 

E = E so1 [GeV] « 0.325 (135) 

for underground experiments in the earth’s mantle 
(p[gm/cm ] = po ~ 3). We find, similarly, that for 
neutrinos of energy E, the solar resonance occurs in 
matter of density 

/9 soi [gm/cm 3 ] sa E[GeV] -1 , (136) 

taking a = 0.0317 from the SNM. 

Likewise, the atmospheric resonance is found for 
neutrinos of energy 

E = E atm [ GeV] w 9.79 , (137) 

also for underground experiments in the earth’s 
mantle. At energy E, it occurs at a density of 

p atm [grn/crn 3 ] sa 29.4/E[GeV] , (138) 

taking cos20 13 = 0.955 from the SNM. 

Because of the close correlation between the 
branch points and the resonances, there is also 
a close correlation between resonances and viable 
approximation schemes. Note, however, that the 
branch point in the a expansion near the solar reso¬ 
nance affects both neutrino and anti-neutrino scat¬ 
tering since its actual location is a = 0. 


D. Discussion 

We have seen in this section that expansions of 
the eigenvalues in a small parameter £' of the SNM 
(where stands for a or sin 2 0i 3 ) must be made 
carefully, since the EV are not analytic everywhere. 

Nevertheless, as noted, first-order Taylor series 
representations of the EV are commonly used to 
simplify the theory. Because these expansions do 
not give an accurate representation the EV for all 
values of the interaction strength A, it is important 
to identify the regions where the theoretical errors 
of the expanded EV might be acceptable and lead 
to accurate representations. We address this in the 
next section. 


VI. APPROXIMATING V{vi v } ) WITH 
EXPANDED EV 

In this section, we begin our assessment of com¬ 
mon procedures used to simplify calculation of the 
oscillation probability by expanding it in one of the 
small parameters of the SNM. Freund observed jij 
that the a-expansion, although useful, could not 
be used near the solar resonance where A = a. 
However, no understanding of the limitations of the 
sin 2 0 13 -expansion appears in the literature. 

We have shown above that the applicability of 
both expansions is limited by the presence of branch 
points in the analytical structure of the eigenvalues. 
The branch point responsible for the failure of the 
a-expansion is located at A = 0 when a = 0, and the 
branch point responsible for the failure of the sin 2 0 33 
expansion is located at A = Aq when sin 2 0 33 = 0, 
where A 0 is defined in Eq. (11291) . 

We make our assessment numerically, com¬ 
paring the oscillation probability calculated from 
Eqs. (I69I80D using the exact eigenvalues to that 
calculated from Eqs. (1691801) using eigenvalues ex¬ 
panded in one of the small parameters of the SNM. 


A. Assessing Oscillation Probabilities 
Expressed in terms of ^-expanded Eigenvalues 

Theoretical errors characterizing an approxima¬ 
tion scheme may emerge numerically only from an 
examination of the dependence of P eM ( Ax, , A) on Ax 
and A. In this section we discuss how we will do this. 

It is convenient to discuss the oscillation pattern 
in terms of the location of the first maxima of the 
functions Ajjo(A^) appearing in the expressions for 
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the partial oscillation probabilities. These peaks oc¬ 
cur at Al = A^\ with 

A { [> = —^— (139) 

2AE[£] 

closely related to the period Pt of jo(Af), 

Pt = 4A^ = — . (140) 

L a m 

Because the exact eigenvalues never cross, the or¬ 
dering of E( is the same as it is in the vacuum, 
namely E 3 > E 2 > E\. It can inferred from this 
that all A E[£] remain positive, and, in addition, that 

AE[2] > AP[3] > 0 

AE[2] > AP[1] > 0 . (141) 

We see, in general, that for small A, AP[1] > AP[3], 
and, for large A, AE[2] > AP[3]. 

It is also clear that AE[2] is always the largest 
eigenvalue difference. Consequently, P 2 is always the 
smallest of the three periods, thus characterizing the 
most rapidly varying Bessel function. The relative 
sizes of Pt are easily worked out in specific cases. In 
the vacuum, 

Pi = 2tt/(1 - a) 

P 2 = 2 t r 

P 3 = 2 tt /a , (142) 

evaluated from differences of the vacuum eigenvalues 
appearing in Eq. (l30l) . 

In the SNM, we find that 

P 3 > Pi for A < A 2 

Pi > P 3 for A > A 2 , (143) 

where 

i 2 = 0.538 (144) 

is the value of A at which P 3 = Pi (AP[ 1 ] = AP[3]). 

With AE[2] the largest eigenvalue difference, P 2 
is always the smallest of the three periods. Thus, 
the value of Al at the first peak of AlJo(A 2 ) is a 
natural scale. 

B. Regions of Maximum Sensitivity 

Sensitivity to the approximation scheme should be 
most manifest within regions where all three Bessel 


functions of P ab (Az J1 A) are of similar size and in¬ 
terfere. This will happen once Al becomes com¬ 
parable to the first peak of its most slowly vary- 

/>. /q\ 

ing A L j 0 (At), which occurs at either A L = A y L ’ or 

A l = A^\ In general, the sensitivity to approx¬ 
imations increases as the distance Al to the most 
distant peak increases. 

In the vacuum, the first peak of the most slowly 
varying AlJo(A^) is always at A^ = A^ 3 \ read¬ 
ily established by the vacuum eigenvalue differ¬ 
ences given in Eq. The corresponding base¬ 

line is obtained from Eqs. (13411391) . l/ 3 )[Km] « 
51707rP[GeV], 

In the SNM, the most slowly varying Bessel func¬ 
tion, established from Eqs. (1143I144|) . is AlJo(A 3 ) 
when A < A 2 and Al. 7 o(A;i ) for A > A 2 . 

1. Regions of Maximum Sensitivity for Fixed A 

Consider first the variation of p( ab ) (Al , A) with 
Al for a given value of A in the SNM. 

According to Eq. (11431) . for A > A 2 , A/jo(A-i) is 
the most slowly varying (having the larger period), 
and the region of maximum sensitivity is 

Al > A< 1} = -4— (145) 

2AP[1] 

with AP[ 1 ] evaluated at A. 

Similarly, for A < A 2 , Alj'o(A 3 ) is the most 
slowly varying. According to Eq. (11391) . its peak 

occurs where Al = 7 r/( 2 A£’[ 3 ]). Thus, the region 
of maximum sensitivity is 

Al > A® = -4— . (146) 

2AP[3] 

with AP[3] again evaluated at A. 

2. Regions of Maximum Sensitivity for Fixed Al 

Consider next the variation of P^ ab \A a L , A) with 
A for given A° in the SNM. 

According to Eq. (11431) . for A > A 2 , Alj'o(Ai) 
is the most slowly varying. According to Eq. (11391) . 
its peak occurs where A 3 = 7 t/( 2 AP[ 1 ]). Thus, the 
region of maximum sensitivity is 

A%] > ^ (147) 
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with AE[ 1 ] evaluated at A. 

Similarly, for A < A 2 , we find from Eq. (11451) that 
Aijo(A 3 ) is the most slowly varying, and the region 
of maximum sensitivity is 

A %] > ^0 ■ ( 148 ) 


VII. NUMERICAL STUDY OF V(v e ->■ z/ M ) 
WITH EXPANDED EV 

Our main interest in the present section is to map 
out the regions where each of the small-parameter 
expansions is capable of simplifying P“ 6 (Al,A). 
We do this by comparing three calculations using 
V{y e —> Vp) taken from our exact Hamiltonian for¬ 
mulation. 

One of these is a calculation of the exact os¬ 
cillation probability obtained in our Hamiltonian 
formulation. For this we use the expressions in 
Eqs. (168I69D evaluated with the exact EV. The other 
two are calculations of our ((-expanded oscillation 
probability for each of the small parameters of the 
SNM. For these we evaluate Eqs. (1681691) using the £- 
expanded EV. For £ = a, we use the a-expanded EV 
given in Eq. (11331) . and for £ = sin 2 $13 we use the 
sin 2 0 i 3 -expanded EV given in Eq. (11311) . The cal¬ 
culation with the ((-expanded EV would, of course, 
coincide with the exact calculation in the vacuum. 
Differences therefore reflect medium effects. 

The extent to which our oscillation probability 
evaluated with one of the ((-expanded EV agrees with 
the exact result indicates regions in which it may 
be possible to obtain, at least in principle, a simple 
((-expanded expression for the oscillation probabil¬ 
ity in good agreement with the exact result. In a 
subsequent paper we make a similar compari¬ 
son between the exact oscillation probability and the 
approximate ones given in Refs. [3, Q. 

From the numerical results we obtain in Ref. H3 
and the present paper, we will be able to identify 
regions in which both ( 1 ) the exact oscillation prob¬ 
ability and the results of Refs. as are in poor agree¬ 
ment; and, ( 2 ) the exact oscillation probability and 
the ((-expanded result are in excellent agreement. 
The regions where both of these conditions are sat¬ 
isfied indicate where it might be possible to improve 
the results found in Refs^jl, S using our Hamilto¬ 
nian formulation. We explore this possibility in yet 
another paper [Is! ]. 

Equations (13411181) provide a means to extrapo¬ 
late the results in any figure to a variety of baseline 


values, medium properties, and neutrino energies 
within the particular regions shown in that figure. 
In the earth’s mantel, where the average Z/N = 1/2, 
Eqs. (11181341) become 

a A 

LlKml = 1.08 x 10 4 -=- 

p[gm/cm 3 ] 

£[GeV] = 30.8 A . (149) 

Plgm/cm J 


A. Ai-Dependence of V(v e — > u fl ) 

We begin our exploration of the extent to which a 
particular ((-expansion is capable of simplifying the 
oscillation probability by examining p( e v) (A l,A) vs 
A l for particular values of A. One value of A is 
chosen near the solar resonance and another near 
the atmospheric resonance. For each choice of £ and 
A, we compare the exact result to Eq. (11301) . 

We first examine V{y e —> v^) below the solar 
resonance, at A = 0 . 0102 . For this value of A, 
AS[3] ~ 0.0294, Eq. (11461) specifies that the ap¬ 
proximate oscillation probability becomes sensitive 

(3) 

to approximations for A^ > A£' « 53. 

Our calculations for A = 0.0102 are shown in 
Fig. [TJ We see from this figure that the a-expanded 
oscillation probability begins to departs from the ex¬ 
act result at large Az, meaning, as expected, that 
the a expansion breaks down in the vicinity of the 
solar resonance. The sensitivity to medium effects 
shows up already for Az ss 20 , which is smaller 
than A^ « 53 estimated using Eq. (11461) . On the 
other hand, Eq. (H30| evaluated with eigenvalues ex¬ 
panded to first order in sin 2 $13 agree well with the 
exact result at large Az showing that the sin 2 $ 13 - 
expanded oscillation probability is capable of pro¬ 
viding an excellent approximation in the vicinity of 
the solar resonance. 

For this small value of A, we find that position 
of the first peak of the exact oscillation probability, 
A l ~ 1.58, coincides almost exactly with the loca¬ 
tion of the peak of the most rapidly varying Bessel 
function, j'o(A 2 ). From Eq. (11491) . we note that the 
oscillation probability at Az = 20 in Fig. [[) where 
the approximate calculation begins to break down in 
the solar resonance region with the a-expanded EV, 
would correspond to a measurement of 105 MeV neu¬ 
trinos propagating in the earth’s mantel at a baseline 
734 Km. 

We next examine V(v e —>■ for A = 0 . 8 , a value 
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FIG. 1 . P £M (Al,A) in the solar-resonance region 
(A = 0.0102) over the interval 15 < Al < 55 for neu¬ 
trinos in matter. Parameters are taken from the SNM. 
Exact result (solid curve). Our expression for the os¬ 
cillation probability evaluated with the a-expanded EV 
(medium-dashed curve). Our expression for the oscil¬ 
lation probability evaluated with the sin 2 #i3-expanded 
EV (long-dashed curve) . 

of A near the atmospheric resonance. Taking the 
exact eigenvalue difference A£7[l] ss 0.328 at A = 
0.8, Eq. (11451) specifies that the desired sensitivity of 
the oscillation probability to approximations should 
become apparent at Al ~ 4.8. Our calculations for 
A = 0.8 are shown in Fig. [ 2 ] 

The exact result shown in Fig. [2] (solid curve) be¬ 
gins to differ from the sin 0 i 3 -expanded result (long- 
dashed curve) at Al ~ 2, which occurs a bit before 
Al « 4.8, where all three Bessel functions fully con¬ 
tribute. Because the long-dashed curve begins to de¬ 
part from the solid curve at large Al, these results 
confirm that P 6 A 1 (Al,A) evaluated with sin 2 613 - 
expanded EV breaks down near the atmospheric res¬ 
onance with medium effects included. The failure of 
the sin (9i 3 -expansion becomes increasingly apparent 
as Al increases to larger Al- 

The fact that the exact result (solid curve) and 
a-expanded result (medium-dashed curve) seem to 
completely overlap demonstrates that P 6 M (Al,A) 
evaluated with the a-expanded EV is capable of be¬ 
coming a completely acceptable approximation near 
the atmospheric resonance. 

For this larger value of A, the first peak of the os¬ 
cillation probability, at Al ~ 1.58, nearly coincides 
with the peak of the most rapidly varying Bessel 
function, j 0 ( A 2 ), at Al ~ 1.5. Applying Eq. (11491) . 
we see that the oscillation probability at the value 
of Al where medium effects begin to become ap¬ 
parent in Fig. [2] would correspond to neutrinos of 
energy EfGeV] a: 4.5 propagating in matter of den¬ 
sity similar to the average density of the entire earth, 


pe/j 



FIG. 2. P e)1 (Al, A) in the atmospheric resonance re¬ 
gion (A = 0 .8) over the interval 0 < Al < 35 for neu¬ 
trinos in matter. Parameters are taken from the SNM. 
Exact result (solid curve). Our expression for the os¬ 
cillation probability evaluated with the a-expanded EV 
(medium-dashed curve). 

p[gm/cm 3 ] = 5.52 at a baseline of 7830 Km (for com¬ 
parison, the average earth radius is 6370 Km). 


B. A-Dependence of V(v e —> 

We next compare oscillation probabilities over 
various ranges of A. For a given range of A , whether 
or not the three Bessel functions maximally interfere 
depends on the choice of Al, which is determined 
by Eqs. (I147I148|) depending on whether A > A 2 or 
A < A 2 , respectively. 

1 . 0 < A < 0.2 

Numerical studies using Eqs. (114711481) show that 
for 0 < A < 0.2 taking Al = 60 is sufficient to 
ensure that the three Bessel functions maximally in¬ 
terfere. In Fig. [3l we compare the oscillation proba¬ 
bilities in this region. These results confirm that the 
expansion in sin 2 $13 is a reasonably good approx¬ 
imation within the solar resonance region, whereas 
the expansion in a is evidently not. 

2. 0.2 < A < A 2 

Numerical studies using Eqs. (I147I148D show that 
for 0 < A < A 2 taking Al = 10 is sufficient to 
ensure that the three Bessel functions maximally in¬ 
terfere. In Fig. [4j we compare the oscillation prob¬ 
abilities in this region. These results confirm that 
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FIG. 3. P bi1 (Al,A) for Al = 60 over the interval 
0 < A < 0.2, for neutrinos in matter. Parameters are 
taken from the SNM. Exact result (solid curve). Our ex¬ 
pression for the oscillation probability evaluated with the 
a-expanded EV (medium-dashed curve). Our expression 
for the oscillation probability evaluated with the sin 2 (9i3- 
expanded EV (long-dashed curve) . 
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FIG. 4. P el1 (Al, A) for An = 10 over the interval 
0.2 < A < A 2 , for neutrinos in matter. Parameters are 
taken from the SNM. Exact result (solid curve). Our ex¬ 
pression for the oscillation probability evaluated with the 
a-expanded EV (medium-dashed curve). Our expression 
for the oscillation probability evaluated with the sin 2 # 13 - 
expanded EV (long-dashed curve) . 

both the expansion in sin 2 $13 and the expansion in 
a are reasonably good approximations here. 

Results shown here apply over the same range of 
neutrino energy, baselines, and medium properties 
as those given in Eqs. (I149D . The extrapolation ap¬ 
plies, of course, only within the region 0.2 < A < A 2 . 

3. A 2 <A< 0.8 

Numerical studies using Eqs. (114711481) show that 
for A 2 < A < 0.8 taking A l = 4 is sufficient to 
ensure that the three Bessel functions maximally in¬ 
terfere. In Fig. [5j we compare the oscillation proba¬ 
bilities over this region. These results confirm that 
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FIG. 5. P E,J (Ai, A) for Al =4 over the interval 
A 2 < A < 0.8, for neutrinos in matter. Parameters 
are taken from the SNM. Exact result (solid curve). Our 
expression for the oscillation probability evaluated with 
the a-expanded EV (medium-dashed curve). Our ex¬ 
pression for the oscillation probability evaluated with the 
sin 2 0i3-expanded EV (long-dashed curve) . 

the expansion in a is a reasonably good approxima¬ 
tion below the atmospheric resonance. The onset 
of the failure of the sin 2 (^-expanded EV near the 
atmospheric resonance begins to become visible for 
A > A2 . Results shown here apply over the same 
range of neutrino energy, baselines, and medium 
properties as those given in Eqs. (I149D . The extrap¬ 
olation applies, of course, only for A 2 < A < 0 . 8 . 

4. 6.8 < A < 1.2 

Numerical studies using Eqs. (114711481) show that 
for 6.8 < A < 1.2 taking Al = 4 is sufficient to 
ensure that the three Bessel functions maximally in¬ 
terfere. In Fig. 0 we compare the oscillation proba¬ 
bilities over this region. These results confirm that 
the expansion in a is a reasonably good approxima¬ 
tion across the atmospheric resonance. The sin 2 # 13 - 
expanded oscillation probability is not shown be¬ 
cause it fails here. Results shown here apply over 
the same range of neutrino energy, baselines, and 
medium properties as those given in Eqs. (I149D . The 
extrapolation applies, of course, only for 0.8 < A < 
1 . 2 . 


5. 1.2 < A < 2.5 

Numerical studies using Eqs. (11471) show that for 
1.2 < A < 2.5 taking A l = 4 is sufficient to en¬ 
sure that the three Bessel functions maximally in¬ 
terfere. In Fig. [71 we compare the oscillation prob- 
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FIG. 6. P eM (Ai,A) for A l = 4 over the interval 
0.8 < A < 1.2, for neutrinos in matter. Parameters 
are taken from the SNM. Exact result (solid curve). Our 
expression for the oscillation probability evaluated with 
the a-expanded EV (medium-dashed curve). Our ex¬ 
pression for the oscillation probability evaluated with the 
sin 2 #i 3 -expanded EV (long-dashed curve) . 





FIG. 7. P bij '(Al, A) vs A for neutrinos in matter over 
the interval 1.2 < A < 2.5 taking Al = 4. Parame¬ 
ters are taken from the SNM. Exact result (solid curve). 
Our expression for the oscillation probability evaluated 
with the a-expanded EV (medium-dashed curve). Our 
expression for the oscillation probability evaluated with 
the sin 2 #i 3 -expanded EV (long-dashed curve) . 


abilities over this region. These results show that 
P 6 M (Al,A) evaluated with the a-expanded EV is 
relatively accurate here. The sin 2 ^-expanded os¬ 
cillation probability agrees with the exact result rea¬ 
sonably well for A > 1.6. 

From Eq. (11491) . we see that with the value = 
4, and in matter of mean density similar to that of 
the average density of the entire earth, p[gm/cm 3 ] = 
5.52, the breakdown of the sin 2 0i 3 -expanded EV 
starts to becomes less visible for A > 1.2, corre¬ 
sponding to neutrinos of energy F[GeV] ~ 6.2 and 
baselines greater than L[Km] = 23, 500 (consider¬ 
ably larger than the diameter of the earth). 



FIG. 8. P £ ^(Al,A) for A l = 35 over the interval 

—0.5 < A < 0.5, for neutrinos and anti-neutrinos in mat¬ 
ter. Parameters are taken from the SNM. Exact result 
(solid curve). Our expression for the oscillation probabil¬ 
ity evaluated with the a-expanded EV (medium-dashed 
curve). Our expression for the oscillation probability 
evaluated with the sin 2 #i 3 -expanded EV (long-dashed 
curve) . 

6. -0.5 < A < 0.5 

In Fig. [BJ we compare the oscillation probability 
calculated with the exact eigenvalues to the oscilla¬ 
tion probability calculated with the expanded eigen¬ 
values for —0.5 < A < 0.5, taking A l = 35. Figure[8] 
confirms the earlier observations that the expansion 
in sin 2 O 13 is a valid approximation within the solar 
resonance region, whereas the expansion in a is not 
for both neutrinos and anti-neutrinos. 


C. Intrinsic Limitations 

It is important to note that when using expanded 
EV, the range of validity of the oscillation proba¬ 
bility P ab (A/,, A) is limited even when it is eval¬ 
uated using the convergent eigenvalue expansion. 
This arises because the error in the phase of the 

trigonometric functions sin(ALAP [i]) grows with 
A^. Eventually, with increasing A/,, the error in 
the phase will lead to an unacceptable error in 

P ab (A L , A). 

When spurious effects of this type would show 
up from a comparison of the exact partial oscilla¬ 
tion probability evaluated with the exact E(_ to the 
exact partial oscillation probability evaluated with 
the convergent expansion of E(. A divergence be¬ 
tween the calculated results of these two calculations 
herald the limit of A^ beyond which the use of ex¬ 
panded eigenvalues breaks down for P aft (Ai, A). 

It is clear from Fig. |T| that the region of validity 
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for both the sin 2 ^-expansion and the a-expansion 
extend out as far as = 35. 


D. Discussion 

The results of this section suggest a natural di¬ 
vision of the full range of A into regions. One sig¬ 
nificant region, which we call the solar resonance 
region, covers the range 0 < A < 0 . 2 . The solar- 
resonance region contains, of course, the solar reso¬ 
nance at A = a. Another is the atmospheric reso¬ 
nance region containing the atmospheric resonance 
at A » cos 2013. This region covers the interval 
A 2 < A < 2.0. 

The region between the solar and atmospheric re¬ 
gions, 0.2 < A < A 2 , is a transition region in the 
sense that within it the expansion in a is improving 
rapidly and the expansion in sin 2 $13 rapidly dete¬ 
riorating. Within this region, the EV expansions 
in sin 2 $13 and a are of comparable accuracy. The 
asymptotic region covers the interval A > 2 . 0 , where 
the EV’s are approaching their asymptotic behavior. 

In future work, we will make a comparison of the 
our oscillation probability evaluated with the exact 
eigenvalues to the approximate oscillation probabil¬ 
ity given in Refs. [3, Q. By comparing those re¬ 
sults to our oscillation probability evaluated with 
^-expanded EV found in the present paper, we will 
be able to identify regions where our Hamiltonian 
formulation might lead to more effective approxi¬ 
mations. Based on this information, we will sub¬ 
sequently present any new approximate results we 
believe be helpful. 


VIII. SUMMARY AND CONCLUSIONS 

We have presented exact, closed-form expressions 
for the neutrino oscillation probabilities in matter 
using our Hamiltonian formulation within the frame¬ 
work of the Standard Neutrino Model assuming 3 
Dirac Neutrinos. Our goal is to benchmark approx¬ 
imate formulations having known difficulties arising 
from expansions commonly used to model neutrino 
and anti-neutrino experiments envisioned for future 
neutrino facilities. 

We have shown explicitly that for small a and 
sin 0 i 3 there are branch points in the analytic struc¬ 
ture of the eigenvalues that lead to singular behavior 
of expansions near the solar and atmospheric reso¬ 
nance. The numerical calculations presented indi¬ 
cate regions in which the small-parameter expan¬ 


sions are guaranteed fail and should be avoided in 
practice. 

Appendix A: Time-evolution operator S(T) in 
Our Hamiltonian Formulation 

Our Hamiltonian formulation was presented in 
Sect. MI Al In this appendix, and in accord with the 
discussion in Sect. IIII Al we give explicit analytic 
formulae for the time-evolution operator S(L) = 

e -xH„A L j n terms of E( and the parameters defining 

H v . 

1. Expression for S(T) in terms of F[t b 

Introducing dimensionless variables, the time- 
evolution operator S(T ) in Eq. (1541) may be written, 

S ab (L) = J2 F e b cos2E e A L 

e 

- i Y Ff sin 2£^A L . (Al) 

1 

Here F^ b , defined in Eq. (f55|) . may be found by eval¬ 
uating Eq. ED at L = 0 and obeys the normaliza¬ 
tion condition 

5 ab = Y F ? b - (A2) 

t 

Using the identity cos 2/3 = 1 — 2 sin 2 /3, we write 
S ab (L) = Y F t b ( l - 2 sin2 E t A L ) 

t 

-iY F ? b sin2^A L 

t 

= Y F ? b ~ 2 F ? b sin2 F(Al 

t 1 

- i Y F i b sin 2F e - a l ■ ( A3 ) 

i 

Now, using the normalization, Eq. (IA2I) . 

S ab {L) = S ab - 2Y F t sin 2 E t A L 
e 

- i Y F i b sin2£ £ A L . (A4) 

l 

Thus, for calculating P ab (L) in Eq. (IMjl . 

Re[S ab (t', t)} =S ab - 2Y M F e b ] sin 2 E t A L 

e 
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+ ^ Im[Fp b ] sin2i5^Ai 

l 

Im[S ab (t',t)} = -2 J2 Irn\ F i b ] sin 2 Eg A L 

l 

- Re i F e b } sm2E e A L . (A5) 

i 

2. Analytic Expressions for Wo[f], W s i„[£], and 
W C oa[^l * n terms °f H v 


+ is\n8c V W sin [i} . (A9) 

with Wq[(\, W cos [i] and W S i n [£] real and indepen¬ 
dent of S cp . Additionally, W s i n \£\ anti-symmetric, 

whereas W CO s[(] and Wq[£] are symmetric, under 
exchange of initial and final states. It follows that 
TrW cos [(-] = TrW S i n [£\ = 0, and therefore 

TrW{£) = TrW 0 [£\ . (A10) 

Then, from Eq. (IA10I) and the definition of Fp b in 
Eq. (l55l) . we find the trace relationship 


In this appendix we give exact, analytic expres¬ 
sions the S(T) in terms of Eg and the parameters 
defining H u . We begin with Eq. (15B1) . which ex¬ 
presses Ff 3 in terms of an operator W[£] and most 
easily accomplish our goal by making the following 
rearrangement of terms in W[£], 

W[£] = UH 0v Ui - UH 0v U f T,{£] 

+ UH 0v U^V + VUH 0v U^ 

+ V 2 -V T,[£] + 1 H[£\ . (A6) 

Here, £[£] is the sum, and n[f] the product, of two 
EV as defined in Eq. (1431441) . 

It is useful to recall that £[£] and II [£], defined in 
this fashion, depend on £ entirely through Eg. In 
particular [l9|, from Eq. (l45l) . 

Z[£\ = -a-E e , (A7) 

and, from Eq. (l48l) . 

n[f] =b + aEg + Eg , (A8) 

with a, b given in Eq. (ED- 

Explicit expressions for W[£] are easily found in 
terms of H „ using Eq. (IA6D . from which it follows 
that its entire dependence on S cp occurs through 
three operators, Ho[f], W cos \£ ] and W S i n [£], 


TrF, = 

m 


= 1 . 


(All) 


a. Analytic Expressions for ITo[f] 

For flavor changing transitions, Vi —> Vj, the ma¬ 
trix elements of W q[£] are found to have the struc¬ 
ture 

< M{fi)\W 0 [£]\M(e) > = W^{£] = c (12) [f?] cos0 23 

< M(t)\W o[£\\M(e) > = W^[£\ = -c^[£] sin0 23 

23 

< M(T)\W 0 [£}\M(fj,) > = Wo [£} 

= c (23) [£} sin 2023 • (A12) 

We use the convention that a quantity 0^ ab ' 1 written 
with parentheses enclosing ab is a number, not a ma¬ 
trix element. As noted, O ab written without paren¬ 
theses surrounding ab is a matrix element. This dis¬ 
tinction may be obvious in the present context, but 
later on this distinction may not be so obvious. 

i ep a er 

Note that W 0 [t] and W 0 [£} are not indepen¬ 
dent. The corresponding coefficients ef' n ' > [£} of for 

-2 

flavor preserving transitions are defined with Ep in 
Eq. (1A8D separated out, 

< M{i)\Wo[£]\M(i) > = Wq[£] 

= E t +c {ii) [£] . (A13) 


W[£\ = W 0 [£] + cosS cp W cos [£} 

I 


The matrix for ITo[f] is thus 


2 / C^[£] C*- 12 ) [£\ COS 023 — c^ 12 ' , [£] sin 023 \ 

W 0 [t\ = lxE l +\ c^[£] cos 0 23 c^[£] c( 23 ) [f] sin 2(9 23 

\ -c (12 )[^] sin 023 c^ 23 ^ [£] sin 2023 c^ 33 ^[£] / 


(A14) 
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Writing the dependence of c^> \(f\ on Eg explicitly, 

c W) M = c « ) +c^ ) E e , 


m „ 

and expressing the results in terms of the combinations of mixing angles 

7j ± ' ) = cos 2 0 i 2 cos 2 023 ± sin 2 0 2 3 sin 2 0 i 2 sin 2 0i 3 


(A15) 


Cl 

C 2 


'~ ±S> = cos 2 0 i 2 ± sin 2 0 i 2 sin 2 0i 3 , 


where C^ was defined earlier in defined in Eq 
changing transitions, 

J12) _ J12) _ « n air, 

Cq — —— — — COS t /13 Sin Zi7i2 


(A16) 

. 63, we find the following exact results. For the flavor- 


and, for the flavor-preserving transitions 


; (23) _ _^ sin 2 Q 12 _ cog 2 g i2 g j n 2 g^ _ ^.( cos 2 g ±3 _ a (j( 1) 

: 1 23) = ^ (cos 2 0 13 - aC ( 2 ~ ] ) ; 


4 n) = a cos 2 0 i 2 cos 2 0 33 
4 n) = — cos 2 0 33 — aC < 2 + ' ) 


(A17) 


(A18) 


•q 22) = a(l — sin 2 0 23 cos 2 0 33 — C^) + A(cos 2 0i 3 cos 2 02 3 

; ( 2 2) _ _|_ cog 2 0 13 g j n 2 g _ £f(+)) _ ^4 _ 


e 23 + a(ci +) -C[ +) )) 


.. - . , , (A19) 

The coefficients c^ 33 ^ follow from c^ 22 ' > by the exchange sin02 3 •<-»• cos02 3 . Consequently, there are only four 
independent matrix elements for Wo[C], two diagonal and two off-diagonal. With the eigenvalues independent 
of 02 3 , it follows that 

~ TT ~ fA(A 

[i\ = Wo ^]|sine 23 ocose 2 3 • (A20) 


b. Analytic Expressions for W s in[t] andW C o S [(] 


The matrices for W S in [P\ and W cos [C\ have the following simple structure 


and 


W sin \i] = a( 12 )[ 


0 — aC 2 )[£] sin 0 23 — aC 2 )[£] cos0 23 


i.( 12 )f 


] sin 0 23 
I cos 023 


0 

-«( 23 )[ 


7 .( 23 ) r 

0 




0 a,( 12 ) 

W cos [£} = ( aC 2 )[£] sin02 3 -a( 23 H 


[P\ sin 0 : 


23 


7.(12) 


[£\ COS 0; 


23 


] sin 2023 — a ( ' 23 ' 1 [(] cos 20 2 3 


aC 2 ) [£] cos 023 —a^ 23 ^[£] cos 2023 a^ 23 ^] sin 202 3 


(A21) 


(A22) 


where 


7 (12) m = ^12) 


,( 12 ) 


a v “'|t| = Ee 


(A23) 


,(23) 


0 = 


(23) , (23) 


l 0 


a 


E e 
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We find the exact results 

= —sin 2^12 sin 0 13 — — sin 20 12 sin 0 13 
a 3 23 ^ = ^ sin 2012 sin 0i 3 
a[, 12) = — ^ cos 2 0 i 2 sin20i 3 
a i 12 ^ = \ sin 20 i 3 — ^ sin 2 0 32 sin 20 i3 . (A24) 

c. Interconnections among Wo[£], W a in{£], and 

w cos [e\ 


Significant simplifications arise from interconnec¬ 
tions among Wo[f], W cos [£], and W S i n [£]. One of 
these follows from Eos. (lA21IA22l) . which show that 

the matrix elements of W cos [I] and W S i n [I] are pro¬ 
portional, namely 

< M(a)\W sin [£}\M(b) > 

= f ab < M (a)\W cos [£)\M (b) > , (A25) 

with the constants of proportionality f ab indepen¬ 
dent of f, 

(° ~ 1 - 1 \ 

fab = 1 0 - cos- 1 2023 . (A26) 

\ 1 cos- 1 2023 0 ) 

Another follows from Eqs. (IA17IA24I) . from which it 
follows that 

a a2) c (12)_ a (12) c (12) =2A , 

a^c^-a^c^ = -K. (A27) 

The quantity K is 

_ a(l - a) 

8 

x cos 0i 3 sin 2012 sin 20 i 3 . (A28) 


3. Expressions for S(T) in terms of S cp 

The dependence of S(T) on the CP violat¬ 
ing phase S cp is very simple and follows from 

Eqs. (1A5IA9I) . noting that Ff 1 = W [£\/D\£\, 
Eq. (1571) . We thus find, 

^ ab 

Re[S ab (t ', f)] = - 2 V -^2-E sin 2 E e A L 

t m 


* ab 

— 2 COS S cp Wcos g j n 2 

t m 

ab 

+ sin 6 cp V w ; in ^ sin2i^AfcA,29) 

e m 


Im[S ab (t',t)\ = -2sin5 cp y sin 2 E e A L 

_y-WyA sin2 E e A L 

i m 

* ab 

E W \£\ 

f° J sin 2EfAj(A30) 

£ m 

The dependence of S(L) on Ee and the remain¬ 
ing parameters of 7A, is given analytically through 
Eqs. (1A14IA21IA22I) . 


Appendix B: Oscillation Probability V{y a —► v b ) 
in our Hamiltonian formulation 

The neutrino oscillation probabilities are obtained 
directly from our expression for the time-evolution 
operator, Eq. (15H) . In the high-energy limit, the 
oscillation probability in Eq. m, expressed in terms 
of Re[Ff b ,] and Im[Ff b ,], is 

V(v a -^■v b ) = 2'y Im[F$,] sin (2AE U ’A L ) 

t>f 

+ 2 y Re[F#} cos {2AFJ U ,A L ) , (Bl) 

i>v 

where 1 < £ < 3, where Al was defined in Eq. (1341) . 
and AEw in Eq. (l44l) . 

Note that Re[F$l] satisfies a normalization condi¬ 
tion 

yReF$ = 5 ab , (B2) 

w 

found by evaluating Eq. (USD at t = t’ and recog¬ 
nizing that V(y a —> V}f) is a real number. Tak¬ 
ing Eq. (IB2I) into account and using the identity 
1 — cos 2/1=2 sin 2 /3, we find an equivalent expres¬ 
sion 

Fiy a —> v b ) = S ab + 2 'y Im[Ff b ,] 

t>t’ 

x sin (2AEw A],) 
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- 4 Y, ReF$ sin 2 {AE U 'A L ) , (B3) 

i>i> 

that bears a striking similarity to the familiar vac¬ 
uum expression with F° b , playing a role analogous to 
J'lf, (as, for example, in Eq. (1) of Ref. |3|). Equa¬ 
tions (IB1IB3I) make use of the fact that a probability 
is purely real. 


1. Properties of F^ b , 

It follows from Eq. El and the observation that 
the energy of a neutrino or antineutrino in matter is 
independent of S cp that F^ b , is symmetric under the 
simultaneous exchange of a,b and £,£', 

F$ = ■ (B4) 

Using in addition the Hermiticity of W[£], Eq. (16U1) 
we find the reflection symmetry 

TT'ba ipab * 

r tf ~ r W 

*tt = F$* . (B5) 

From this, along with Eq. (IB4I) . it follows that 
ImFffi is antisymmetric under the exchange of ei¬ 
ther (a,b) or (£,£'), 

ImF$ = -ImFfi = -ImFfe , (B6) 

whereas ReFffi, is symmetric, 

ReF$ = ReF e ab = ReF% . (B7) 

There are no symmetries connecting the F£ b , of 
neutrinos to those of antineutrinos because these en¬ 
ergies are, in general, different. This is not the case, 
however, in the vacuum for theories invariant under 
CPT. 


2. General Expressions for Fffl 


The quantity F?}, 


w 


ab 

tt’i 


is most easily obtained from 


wfg, = < M{b)\W[£]\M(a) > 

x < M(b)\W[£']\M(a) >* . (B8) 

Analytic formulae for w^\, are easily obtained in 

terms of the parameters of F£ v using W[t] given in 
Eq. (IA9D . Equation (1731) then gives Fgfi as 


jpab _ w gg 1 

r w ~ ~ -s- 

D[i]D{£'} 


(B9) 


Because W[£] consists of three terms, one propor¬ 
tional to smAcp, one proportional to cos<5 cp , and one 
independent of S cp , the dependence of w ab [£, £'] on 
S cp can be expressed a priori through the five oper¬ 
ators, 

wu' = w o\f + cos 6 cp w“ b sU , 

+ cos 2 S cp w^ sHe , + i{svo8 cp wf inU , 

+ sin S cp cos 8 C pVJ S i n x CO sW > (BIO) 

each of which is uniquely determined by 
Eqs. (IA9IB8I) . Clearly, just as for the Wj {(], 
the matrix elements of u>i [£, £'} are real and 
independent of 8 cp . 

Using the hermiticity of W${£\ and W cos [£], and 
the anti-hermiticity of W S i n [£\, we then find 

, ± ab ± ab ^ ab ± ab 

<U' = m 0 [t] + W sin [£]W sm [l'] 

. * ab ± ab ^ ab ^ ab 

<lw = W Hn [£\W 0 [U] - IU 0 [£\W sin [£'} 

. i ab ± ab ~ ab * ab 

<osU' = W^WWo [£'} + W 0 [£}W cos [£'} 

~ ab ± ab ± ab ± ab 

= W cos m cos [F] - W sin [t]W sin [t'] 

^ ab ^ ab 

<l*cosW = W sin [£]W cos [£'} 

± ab ± ab 

-W cos m s in[t'} ■ (Bll) 

Note that Eqs. (IA25IA26I) require that w S i nxcos van¬ 
ish identically, 

<inxcos = 0 , (B12) 

so this term need not be considered further. Note 
also that the dependence of wf b e , on £ and £' arises 

iL i ab 

entirely from the eigenvalues Eg and Eg> , as W [£} 
depends on £ entirely through Eg. Finally, we will 
find it useful to define w^ b [£} in analogy to W[t] in 
Eq. (l58l) . 

< b [l]=< 6 [3,2] 

wf[2}=wf[3,l} 

< 6 [3]=< 6 [2,1] . (B13) 

These are the only three wf b [£, £'} needed because of 
the restriction £ > £' in Eq. (IB3D . 


3. Analytic expressions for Wsin[£,£'] 

* ab n ab 

From the symmetries of W 0 [f], W cos [t], and 

^ ab 

W sin [(], we see from Eq. (|B 1 1|) that is odd 
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under the exchange of either a, b or t, £'. The term 
w sine i> therefore vanishes for a = b [and for £ = £'], 
but again recall that restrictions on the sums is such 
that wf b , contributes only for i > £'. 

Using the general results in Eqs. (IA14IA21IA22I) . 
we find w g b nU , given in Eq. (IB 11 11 is, as a matrix, 


1 


0 A —A 


w S inW = - sin 26*23 I -A 0 B \ ,(B14) 
A -B 0 


where 


and 


A = -(a^[£]c^[£'} - c^[£]a^[£']) 

= -(a^c^-a^c^)AE w (B15) 


B = 2(a* 23 >[<]c< 23 >[^ - c^[£]a^[£']) 

= 2(4 23) c( 23) -af ) cf 3) )A^«' • (B16) 


We see from Eq. (IA27I) that that A = B, so from 
Eq. (IB 141) we find the following simple expression, 


ab 


wfinW = K Sin 2023 A Ewe': 
and correspondingly 

= Ksm29 23 AE[e}e a J n , 

where e s i n is the anti-symmetric matrix 



(B17) 


(B18) 


(B19) 


4. Analytic expressions for wf b g2 , i and h 


ab 

We also see from the symmetries of W 0 [i\, 

* ab ~ ab 

W cos [£] , and W sin [£] , that w^ w , wf oM , , and w ab os2 u , 
are even under the exchange of a, b and/or £,£'. 
Since the dependence of all wf b ?, on £ and £’ arises 

entirely from the eigenvalues E( and Ee>, it follows 
that the terms symmetric in £ and £' must be func¬ 
tions of symmetric combinations of Ee and E&. 

There are only three irreducible symmetric func¬ 
tions of Et and Ef. a constant, the sum Ee + Ef, 
and the product E(E(i. just as for W[£] in Eq. (IA6I) . 
Thus, we anticipate that the entire dependence of 
Wg b [£], w“ b s [£], and w “ b s2 {£] on £ will occur through 
T,[£] and II [£\ in Eqs. (IA7IA8I) . and hence through 
E t , as for W[£]. 


In the following discussion, it is important to re¬ 
call our convention that a quantity, such as u/ 0,6 - 1 [£], 
written with parentheses enclosing ab is a num¬ 
ber, whereas w ab [£] written without parentheses sur¬ 
rounding ab is the element of a matrix w[£\, 

w ab [£] = < M{b)\w[£]\M(a) > . (B20) 


a. Matrix Elements of w cos 2 

Using Eq. (IB 111) with Eqs. (IA21IA22I) . the matrix 
for w cos 2 is found to be 


= sin 2 29 


23 

0 0 


0 


x | 0 W {22 1 ~W {22 1 

■ COS Z ,COS z 

n (22) (22) 

0 —W 2 W o 

cos z cos z 


, (B21) 


where 


ASne = » <23| [<]o (23) K] 


(23)2 (23) (23) 


x 0 


Aa^a^iEe + Ee,) 


(B22) 


(B23) 


+ a ( 23 ) 2 EeE e ■ 

Correspondingly, 

(22) rn-T (23)2 (23) (23) vr/ 

w coM?]= a o + a o a\ ’ T, [t 

+ 4 23)2 n[£]. 

Note the somewhat subtle notational distinction 
between the operator w cos 2 and the coefficients w 
in terms of which it is defined. The latter is in¬ 
dicated by parentheses that surrounding the super¬ 
scripts. By contrast, superscripts without parenthe¬ 
ses, as in wf > 1 indicate the transition v a —> zxf,, i.e., 
wf b =< M(b)\wi\M(a) >. 


b. Matrix Elements of w c 


We find only three independent, non-vanishing 
matrix elements for w cos , one diagonal element and 
two off-diagonal elements. As a matrix, 


w C os = sin 2023 


/ 

0 

(12) 

Wcos 

(12) 
- Wcos 


,,,(12) 

,„(22) 

■of 23 ') 


Wcos 

Wcos 

Wcos 


"112) 

of 23 ') 

„„(33) 

V 

Wqos 

W C os 

Wcos 


(B24) 


The two independent off-diagonal elements, 

^ = ^(« (12) Mc (12 ) M+c( 12 nfla (12 ^^]) 

= -{a^[£]c^[£'} +c^[£]a^[£']) 
x cos 2023 , (B25) 
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follow immediately from the structure of W 0 [l'\ and 

±. ab __ 

W cos [(] given in Eqs. (IA14IA22I) . respectively. 
Taking w^cos from Eqs. (IA14IA22IB 111) , we obtain 

■»S«. = -« (23> M(4 2 . + <= (22| [C]) 

- (e] + c (22) [£])a (23) [f'] . (B26) 

Similarly, for the diagonal matrix element D , 

+ (i?“ + c (33) M)a (23) M] . (B27) 

The element vanishes as a consequence of 

22 ( 11 ) 

w cos [(] = 0. 

Correspondingly, for the off-diagonal matrix ele¬ 
ments we find, 

»£ 2, M=«r4 12, + 5(» ( . I2> 4 ,2, +a< I2) cS I2 >) 
x EM + a£ 12) c[ 12) nM 

w£ 23) M = -2cos26»23«o 23) c|, 23) 

-co S 20 2 3(a^cf ) +4 23) c^) 

x £[£] - 2cos26» 23 af 3) 4 23) n[£] , (B28) 

which simplify somewhat by using Eqs. (IA27IA28I) . 
For the diagonal matrix elements, 

(22)o (23) (22) / (23) (22) . (23) (22) n 

W co S ’W = ~ 2a 0 c 0 ) ’ + O.Q 4 ; ) 

x £[f] + 2(a{, 23) - af 3) 4 22) ) n[£] 

- af 3) S[^] n[£] - 4 23) £[£] 2 (B29) 

and 

^ 3) M=2a( 23) 4 33) + (4 23 )cr + ar c r) 

x EM - 2(a[) 23) - af 3) cf 3) ) n[£] 

+ af 3) £[£] IIM + 4 23) £[£] 2 . (B30) 

The matrix elements iCcos'* M and wioV M are not in- 
dependent. The reason is that w y cos M may be found 

(oo) 

from Wcos M by making the replacement sin $23 ^ 
cos 0 2 3 [Eq. (IA20I) ] and by flipping the overall sign. 
Recall that under this replacement, c^ 22 -* O c^ 33 - 1 and 
that a ^ 23 ^ is independent of sin0 23 . 

c. Matrix Elements of wo 

We find four independent, non-vanishing matrix 
elements for wq, two diagonal and two off-diagonal 


elements. As a matrix, 


w 0 = 


/ 

(11) 

W K q 

(12) 

^0 

(13) 

Wo 


,,,(12) 

,„( 22 ) 

,,,(23) 


W 0 

Wq 

Wq 

V 

„„(13) 

„„(23) 

„„(33) 

Wq 

w 0 

Wq 


(B31) 


The off-diagonal elements, 


w, 


.(12) 

0«' 

= 

c (!2) 

Mc (12) M] 

cos 2 

023 


+ 


Ma (12) M] 

sin 2 

023 

,(13) 

0«' 

= 

c (!2) 

Mc (12) M] 

sin 2 

023 


+ 

a™ 

M« (12) M] 

cos' 

' 023 

O ' 

to 

^6. W 

= 

c (23) 

[l]cW[t'] 

sin 2 

2023 


+ 



5 



(B32) 

n ab 

follow immediately from the structure of W 0 \t'\ 

~ ab __ 

and W sin \(] given in Eqs. (IA14IA21I) . respec¬ 
tively. The diagonal elements Wq 1 ^ are found from 
Eqs. (1A14IB11D , 


w, 


(nn 


•> — + c 


(nn) r 


X (E t , + c (nn) M]) . 


(B33) 


Correspondingly, for the off-diagonal matrix ele¬ 
ments, 

„,(!2) m __ (12)2 2 „ (12)2 - 2 n 

w 0 [i\ — c 0 cos e 2 3 + a 0 sm y 23 

+ (cq 12) C 3 12 ^ cos 2 0 23 + a[ ) 12 ' ) aj 12 ' > sin 2 0 23 )£[£] 
+ (c[ 12 ^ 2 cos 2 0 23 + a[ 12)2 sin 2 0 23 ) IT[£] 
w o [IJ — Cg sm U-23 + a 0 cos 023 

+ (cq 12) C 3 12) sin 2 0 23 + ao 12 ^Uj 12 ' ) cos 2 0 2 3)£M 
+ (4 12)2 sin 2 023 + a^ 12 '* 2 cos 2 0 2 3) IIM 

w 0 M = c 0 Sm 2023 + % 

+ (4 23) 4 23) +4 23) cf 3) sin 2 20 23 ) SM 
+ (cf 3)2 sin 2 023 + <4 23)2 ) nM . (B34) 

For the diagonal matrix elements, 

4 nn) M = 4” n)2 - (2<> n) - 4 nn)2 ) n [t\ 

+4 nn) 4 nn) sm+ c { r ] m m 

+ nM 2 . (B35) 

As for w cos , the matrix elements Wq 22 ^ [i] and Wg 33 -* [t\ 
are not independent. The same is true for Wq [l\ 
and ; \i\ . In both cases the matrix elements may 
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be found from one another by making the replace¬ 
ment sin 6*23 -H- cos 023 [Eq. (IA20I) ]. Recall that un¬ 
der such a replacement c*- 22 ) F4 and that c^ 12 ^ 
and a^ 23 ) are independent of sin 0 2 3 . 

Writing the dependence of the coefficients 
on Ei explicitly, 


(mn) , (mn) 

= K,o + W U Ee 

. ( mn ) 

+ W\.,2 E i ■ 


(B36) 


5. Analytic expression for Trri[F^l] and Re\F"j ‘,] 


A general form for the real and imaginary parts of 
Fjjfy in terms of the four operators, (w o, w S in , w cos , 
w cos 1 2 3 4 5 6 7 8 9 10 ) is found from Eq. (l73l) and Eq. (IB 101) . For 
Im[Fj$\, using Eq. (IB17I) . we find 


Im[F$] = sin S cp sin f 
D[t]D[F] 

= i sin 5 cp sin 20 23 AT®/;, 


A Ew 

x —- E - . 

D[i]D[t'] 

Then, using the easily-verified result 
A Eu, e«; 

b[i]b[i'\ b 


(B37) 


(B38) 


we find 

ImF$ = - sin S cp e a J n e^' n ^-^ cos0i 3 

8 D 

x sin 2012 sin 2013 sin 2023 • (B39) 

For Re[Ffft,] we find 

Re[F#\ = + cos S cp W *°l u ' 

D[e)D[£'] D[£]D[F] 

w ab 

+ cos 2 S cp cos l 11 ' . (B40) 

D[i]D[i'} 


Equations (157IB38j) show that Ffy plays a role 
for V(v a —> Vb) similar to the one that Fj? b 
plays for S(T). Using Eq. (IB10I) . it follows from 
Eqs. (IB37IB40I) that ImFj? b , is odd, and ReFffl even, 
under S cp —> —5 cp . 
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